Genomic prediction of blood biomarkers of metabolic disorders in Holstein cattle using parametric and nonparametric models.

BACKGROUND: Metabolic disturbances adversely impact productive and reproductive performance of dairy cattle due to changes in endocrine status and immune function, which increase the risk of disease. This may occur in the post-partum phase, but also throughout lactation, with sub-clinical symptoms. Recently, increased attention has been directed towards improved health and resilience in dairy cattle, and genomic selection (GS) could be a helpful tool for selecting animals that are more resilient to metabolic disturbances throughout lactation. Hence, we evaluated the genomic prediction of serum biomarkers levels for metabolic distress in 1353 Holsteins genotyped with the 100K single nucleotide polymorphism (SNP) chip assay. The GS was evaluated using parametric models best linear unbiased prediction (GBLUP), Bayesian B (BayesB), elastic net (ENET), and nonparametric models, gradient boosting machine (GBM) and stacking ensemble (Stack), which combines ENET and GBM approaches. RESULTS: The results show that the Stack approach outperformed other methods with a relative difference (RD), calculated as an increment in prediction accuracy, of approximately 18.0% compared to GBLUP, 12.6% compared to BayesB, 8.7% compared to ENET, and 4.4% compared to GBM. The highest RD in prediction accuracy between other models with respect to GBLUP was observed for haptoglobin (hapto) from 17.7% for BayesB to 41.2% for Stack; for Zn from 9.8% (BayesB) to 29.3% (Stack); for ceruloplasmin (CuCp) from 9.3% (BayesB) to 27.9% (Stack); for ferric reducing antioxidant power (FRAP) from 8.0% (BayesB) to 40.0% (Stack); and for total protein (PROTt) from 5.7% (BayesB) to 22.9% (Stack). Using a subset of top SNPs (1.5k) selected from the GBM approach improved the accuracy for GBLUP from 1.8 to 76.5%. However, for the other models reductions in prediction accuracy of 4.8% for ENET (average of 10 traits), 5.9% for GBM (average of 21 traits), and 6.6% for Stack (average of 16 traits) were observed. CONCLUSIONS: Our results indicate that the Stack approach was more accurate in predicting metabolic disturbances than GBLUP, BayesB, ENET, and GBM and seemed to be competitive for predicting complex phenotypes with various degrees of mode of inheritance, i.e. additive and non-additive effects. Selecting markers based on GBM improved accuracy of GBLUP.

Benchmarking machine learning and parametric methods for genomic prediction of feed efficiency-related traits in Nellore cattle.

Genomic selection (GS) offers a promising opportunity for selecting more efficient animals to use consumed energy for maintenance and growth functions, impacting profitability and environmental sustainability. Here, we compared the prediction accuracy of multi-layer neural network (MLNN) and support vector regression (SVR) against single-trait (STGBLUP), multi-trait genomic best linear unbiased prediction (MTGBLUP), and Bayesian regression (BayesA, BayesB, BayesC, BRR, and BLasso) for feed efficiency (FE) traits. FE-related traits were measured in 1156 Nellore cattle from an experimental breeding program genotyped for ~ 300 K markers after quality control. Prediction accuracy (Acc) was evaluated using a forward validation splitting the dataset based on birth year, considering the phenotypes adjusted for the fixed effects and covariates as pseudo-phenotypes. The MLNN and SVR approaches were trained by randomly splitting the training population into fivefold to select the best hyperparameters. The results show that the machine learning methods (MLNN and SVR) and MTGBLUP outperformed STGBLUP and the Bayesian regression approaches, increasing the Acc by approximately 8.9%, 14.6%, and 13.7% using MLNN, SVR, and MTGBLUP, respectively. Acc for SVR and MTGBLUP were slightly different, ranging from 0.62 to 0.69 and 0.62 to 0.68, respectively, with empirically unbiased for both models (0.97 and 1.09). Our results indicated that SVR and MTGBLUBP approaches were more accurate in predicting FE-related traits than Bayesian regression and STGBLUP and seemed competitive for GS of complex phenotypes with various degrees of inheritance.

Genotype-by-environment interactions for feed efficiency traits in Nellore cattle based on bi-trait reaction norm models.

BACKGROUND: Selecting animals for feed efficiency directly impacts the profitability of the beef cattle industry, which contributes to minimizing the environmental footprint of beef production. Genetic and environmental factors influence animal feed efficiency, leading to phenotypic variability when exposed to different environmental conditions (i.e., temperature and nutritional level). Thus, our aim was to assess potential genotype-by-environment (G × E) interactions for dry matter intake (DMI) and residual feed intake (RFI) in Nellore cattle (Bos taurus indicus) based on bi-trait reaction norm models (RN) and evaluate the genetic association between RFI and DMI across different environmental gradient (EG) levels. For this, we used phenotypic information on 12,958 animals (young bulls and heifers) for DMI and RFI recorded during 158 feed efficiency trials. RESULTS: The heritability estimates for DMI and RFI across EG ranged from 0.26 to 0.54 and from 0.07 to 0.41, respectively. The average genetic correlations (± standard deviation) across EG for DMI and RFI were 0.83 ± 0.19 and 0.81 ± 0.21, respectively, with the lowest genetic correlation estimates observed between extreme EG levels (low vs. high) i.e. 0.22 for RFI and 0.26 for DMI, indicating the presence of G × E interactions. The genetic correlation between RFI and DMI across EG levels decreased as the EG became more favorable and ranged from 0.79 (lowest EG) to 0.52 (highest EG). Based on the estimated breeding values from extreme EG levels (low vs. high), we observed a moderate Spearman correlation of 0.61 (RFI) and 0.55 (DMI) and a selection coincidence of 53.3% and 40.0% for RFI and DMI, respectively. CONCLUSIONS: Our results show evidence of G × E interactions on feed efficiency traits in Nellore cattle, especially in feeding trials with an average daily gain (ADG) that is far from the expected of 1 kg/day, thus increasing reranking of animals.

Integrating on-farm and genomic information improves the predictive ability of milk infrared prediction of blood indicators of metabolic disorders in dairy cows.

BACKGROUND: Blood metabolic profiles can be used to assess metabolic disorders and to evaluate the health status of dairy cows. Given that these analyses are time-consuming, expensive, and stressful for the cows, there has been increased interest in Fourier transform infrared (FTIR) spectroscopy of milk samples as a rapid, cost-effective alternative for predicting metabolic disturbances. The integration of FTIR data with other layers of information such as genomic and on-farm data (days in milk (DIM) and parity) has been proposed to further enhance the predictive ability of statistical methods. Here, we developed a phenotype prediction approach for a panel of blood metabolites based on a combination of milk FTIR data, on-farm data, and genomic information recorded on 1150 Holstein cows, using BayesB and gradient boosting machine (GBM) models, with tenfold, batch-out and herd-out cross-validation (CV) scenarios. RESULTS: The predictive ability of these approaches was measured by the coefficient of determination (R2). The results show that, compared to the model that includes only FTIR data, integration of both on-farm (DIM and parity) and genomic information with FTIR data improves the R2 for blood metabolites across the three CV scenarios, especially with the herd-out CV: R2 values ranged from 5.9 to 17.8% for BayesB, from 8.2 to 16.9% for GBM with the tenfold random CV, from 3.8 to 13.5% for BayesB and from 8.6 to 17.5% for GBM with the batch-out CV, and from 8.4 to 23.0% for BayesB and from 8.1 to 23.8% for GBM with the herd-out CV. Overall, with the model that includes the three sources of data, GBM was more accurate than BayesB with accuracies across the CV scenarios increasing by 7.1% for energy-related metabolites, 10.7% for liver function/hepatic damage, 9.6% for oxidative stress, 6.1% for inflammation/innate immunity, and 11.4% for mineral indicators. CONCLUSIONS: Our results show that, compared to using only milk FTIR data, a model integrating milk FTIR spectra with on-farm and genomic information improves the prediction of blood metabolic traits in Holstein cattle and that GBM is more accurate in predicting blood metabolites than BayesB, especially for the batch-out CV and herd-out CV scenarios.

Heteroscedastic Reaction Norm Models Improve the Assessment of Genotype by Environment Interaction for Growth, Reproductive, and Visual Score Traits in Nellore Cattle.

The assessment of the presence of genotype by environment interaction (GxE) in beef cattle is very important in tropical countries with diverse climatic conditions and production systems. The present study aimed to assess the presence of GxE by using different reaction norm models for eleven traits related to growth, reproduction, and visual score in Nellore cattle. We studied five reaction norm models (RNM), fitting a linear model considering homoscedastic residual variance (RNM_homo), and four models considering heteroskedasticity, being linear (RNM_hete), quadratic (RNM_quad), linear spline (RNM_l-l), and quadratic spline (RNM_q-q). There was the presence of GxE for age at first calving (AFC), scrotal circumference (SC), weaning to yearling weight gain (WYG), and yearling weight (YW). The best models were RNM_l-l for YW and RNM_q-q for AFC, SC, and WYG. The heritability estimates for RNM_l-l ranged from 0.07 to 0.20, 0.42 to 0.61, 0.24 to 0.42, and 0.47 to 0.63 for AFC, SC, WYG, and YW, respectively. The heteroskedasticity in reaction norm models improves the assessment of the presence of GxE for YW, WYG, AFC, and SC. Additionally, the trajectories of reaction norms for these traits seem to be affected by a non-linear component, and selecting robust animals for these traits is an alternative to increase production and reduce environmental sensitivity.

Integrating genome-wide association study and pathway analysis reveals physiological aspects affecting heifer early calving defined at different ages in Nelore cattle.

Heifer early calving (HC) plays a key role in beef cattle herds' economic sustainability and profitability by reducing production costs and generation intervals. However, the genetic basis of HC in Nelore heifers at different ages remains to be well understood. In this study, we aimed to perform a multi-trait weighted single-step genome-wide association (MT w-ssGWAS) to uncover the genetic mechanism involved in HC at 24 (HC24), 26 (HC26), 28 (HC28), and 30 (HC30) months of age in Nelore heifers. The MT w-ssGWAS pointed out four shared windows regions for HC24, HC26, HC28, and HC30 on BTA 5, 6, 14, and 16, explaining a larger proportion of genetic variation from 9.2% for HC30 to 10.6% for HC28. The shared regions harbored candidate genes related with the major gatekeeper for early puberty onset by controlling metabolic aspects related to homeostasis, reproductive, and growth (IGF1, PARPBP, PMCH, GNRHR, LYN, TMEM68, PLAG1, CHCHD7, KISS1, GOLT1A, and PPP1R15B). The MT w-ssGWAS and pathway analysis highlighted differences in physiological processes that support complex interactions between the gonadotropic axes, growth aspects, and sexual precocity in Nelore heifers, providing useful information for genetic improvement and management strategies.

Meta-analysis across Nellore cattle populations identifies common metabolic mechanisms that regulate feed efficiency-related traits.

BACKGROUND: Feed efficiency (FE) related traits play a key role in the economy and sustainability of beef cattle production systems. The accurate knowledge of the physiologic background for FE-related traits can help the development of more efficient selection strategies for them. Hence, multi-trait weighted GWAS (MTwGWAS) and meta-analyze were used to find genomic regions associated with average daily gain (ADG), dry matter intake (DMI), feed conversion ratio (FCR), feed efficiency (FE), and residual feed intake (RFI). The FE-related traits and genomic information belong to two breeding programs that perform the FE test at different ages: post-weaning (1,024 animals IZ population) and post-yearling (918 animals for the QLT population). RESULTS: The meta-analyze MTwGWAS identified 14 genomic regions (-log10(p -value) > 5) regions mapped on BTA 1, 2, 3, 4, 7, 8, 11, 14, 15, 18, 21, and 29. These regions explained a large proportion of the total genetic variance for FE-related traits across-population ranging from 20% (FCR) to 36% (DMI) in the IZ population and from 22% (RFI) to 28% (ADG) in the QLT population. Relevant candidate genes within these regions (LIPE, LPL, IGF1R, IGF1, IGFBP5, IGF2, INS, INSR, LEPR, LEPROT, POMC, NPY, AGRP, TGFB1, GHSR, JAK1, LYN, MOS, PLAG1, CHCD7, LCAT, and PLA2G15) highlighted that the physiological mechanisms related to neuropeptides and the metabolic signals controlling the body's energy balance are responsible for leading to greater feed efficiency. Integrated meta-analysis results and functional pathway enrichment analysis highlighted the major effect of biological functions linked to energy, lipid metabolism, and hormone signaling that mediates the effects of peptide signals in the hypothalamus and whole-body energy homeostasis affecting the genetic control of FE-related traits in Nellore cattle. CONCLUSIONS: Genes and pathways associated with common signals for feed efficiency-related traits provide better knowledge about regions with biological relevance in physiological mechanisms associated with differences in energy metabolism and hypothalamus signaling. These pleiotropic regions would support the selection for feed efficiency-related traits, incorporating and pondering causal variations assigning prior weights in genomic selection approaches.

Real-time milk analysis integrated with stacking ensemble learning as a tool for the daily prediction of cheese-making traits in Holstein cattle.

Cheese-making traits in dairy cattle are important to the dairy industry but are difficult to measure at the individual level because there are limitations on collecting phenotypic information. Mid-infrared spectroscopy has its advantages, but it can only be used during monthly milk recordings. Recently, in-line devices for real-time analysis of milk quality have been developed. The AfiLab recording system (Afimilk) offers significant benefits as phenotypes can be collected from each cow at each milking session. The objective of this study was to assess the potential of integrating AfiLab real-time milk analyzer measures with the stacking ensemble learning technique using heterogeneous base learners for the in-line daily monitoring of cheese-making traits in Holstein cattle with a view to developing a precision livestock farming system for monitoring the technological quality of milk. Data and samples for wet-laboratory analyses were collected from 499 Holstein cows belonging to 2 farms where the AfiLab system was installed. The traits of concern were 9 milk coagulation traits [3 milk coagulation properties (MCP), and 6 curd firming traits (CFt)], and 7 cheese-making traits [3 cheese yield (CY) traits, and 4 milk nutrient recovery in the curd (REC) traits]. The near-infrared AfiLab spectral data and on-farm information (days in milk and parity) were used to assess the predictive ability of different statistical methods [elastic net (EN), gradient boosting machine (GBM), extreme gradient boosting (XGBoost), and artificial neural network (ANN)] across different cross-validation scenarios. These statistical methods were considered the base learners, which were then combined in a stacking ensemble learning. Results indicate that including information on the cows (days in milk and parity) in the AfiLab infrared prediction increased its accuracy by 10.3% for traditional MCP, 13.8% for curd firming, 9.8% for CY, and 11.2% for REC traits compared with those obtained from near-infrared AfiLab alone. The statistical approaches exhibited high prediction accuracies (R2) averaged across the cross-validation scenarios for traditional MCP (0.58 for ANN, 0.55 for EN and GBM, 0.52 for XGBoost, and 0.62 for stacking ensemble), CFt (0.55 for ANN, 0.54 for EN and GBM, 0.53 for XGBoost, and 0.61 for stacking ensemble), and similar R2 averages for CY and REC (0.55 for ANN, 0.54 for EN and GBM, 0.53 for XGBoost, and 0.61 for stacking ensemble). The ANN approach was more accurate than the other base learners (EN, GBM, and XGBoost) and improved accuracy across cross-validation scenarios on average by 7% for traditional MCP, 5% for CFt, 8% for CY, and 7% for REC. The stacking ensemble method improved prediction accuracy by 3% to 31% for traditional MCP, 2% to 26% for CFt, 1% to 38% for CY traits, and 2% to 27% for REC traits compared with the base learners. The prediction accuracies of the different approaches evaluated tended to decrease from the 10-fold cross-validation to the independent validation scenario, although there was a smaller reduction in prediction accuracy with the stacking ensemble learning technique across all the cross-validation scenarios. Our results show that combining in-line on-farm information with stacking ensemble machine learning represents an effective alternative for obtaining robust daily predictions of milk cheese-making traits.

Evaluating the performance of machine learning methods and variable selection methods for predicting difficult-to-measure traits in Holstein dairy cattle using milk infrared spectral data.

Fourier-transform infrared (FTIR) spectroscopy is a powerful high-throughput phenotyping tool for predicting traits that are expensive and difficult to measure in dairy cattle. Calibration equations are often developed using standard methods, such as partial least squares (PLS) regression. Methods that employ penalization, rank-reduction, and variable selection, as well as being able to model the nonlinear relations between phenotype and FTIR, might offer improvements in predictive ability and model robustness. This study aimed to compare the predictive ability of 2 machine learning methods, namely random forest (RF) and gradient boosting machine (GBM), and penalized regression against PLS regression for predicting 3 phenotypes differing in terms of biological meaning and relationships with milk composition (i.e., phenotypes measurable directly and not directly in milk, reflecting different biological processes which can be captured using milk spectra) in Holstein-Friesian cattle under 2 cross-validation scenarios. The data set comprised phenotypic information from 471 Holstein-Friesian cows, and 3 target phenotypes were evaluated: (1) body condition score (BCS), (2) blood ß-hydroxybutyrate (BHB, mmol/L), and (3) κ-casein expressed as a percentage of nitrogen (κ-CN, % N). The data set was split considering 2 cross-validation scenarios: samples-out random in which the population was randomly split into 10-folds (8-folds for training and 1-fold for validation and testing); and herd/date-out in which the population was randomly assigned to training (70% herd), validation (10%), and testing (20% herd) based on the herd and date in which the samples were collected. The random grid search was performed using the training subset for the hyperparameter optimization and the validation set was used for the generalization of prediction error. The trained model was then used to assess the final prediction in the testing subset. The grid search for penalized regression evidenced that the elastic net (EN) was the best regularization with increase in predictive ability of 5%. The performance of PLS (standard model) was compared against 2 machine learning techniques and penalized regression using 2 cross-validation scenarios. Machine learning methods showed a greater predictive ability for BCS (0.63 for GBM and 0.61 for RF), BHB (0.80 for GBM and 0.79 for RF), and κ-CN (0.81 for GBM and 0.80 for RF) in samples-out cross-validation. Considering a herd/date-out cross-validation these values were 0.58 (GBM and RF) for BCS, 0.73 (GBM and RF) for BHB, and 0.77 (GBM and RF) for κ-CN. The GBM model tended to outperform other methods in predictive ability around 4%, 1%, and 7% for EN, RF, and PLS, respectively. The prediction accuracies of the GBM and RF models were similar, and differed statistically from the PLS model in samples-out random cross-validation. Although, machine learning techniques outperformed PLS in herd/date-out cross-validation, no significant differences were observed in terms of predictive ability due to the large standard deviation observed for predictions. Overall, GBM achieved the highest accuracy of FTIR-based prediction of the different phenotypic traits across the cross-validation scenarios. These results indicate that GBM is a promising method for obtaining more accurate FTIR-based predictions for different phenotypes in dairy cattle.

Integrating genomic and infrared spectral data improves the prediction of milk protein composition in dairy cattle.

BACKGROUND: Over the past decade, Fourier transform infrared (FTIR) spectroscopy has been used to predict novel milk protein phenotypes. Genomic data might help predict these phenotypes when integrated with milk FTIR spectra. The objective of this study was to investigate prediction accuracy for milk protein phenotypes when heterogeneous on-farm, genomic, and pedigree data were integrated with the spectra. To this end, we used the records of 966 Italian Brown Swiss cows with milk FTIR spectra, on-farm information, medium-density genetic markers, and pedigree data. True and total whey protein, and five casein, and two whey protein traits were analyzed. Multiple kernel learning constructed from spectral and genomic (pedigree) relationship matrices and multilayer BayesB assigning separate priors for FTIR and markers were benchmarked against a baseline partial least squares (PLS) regression. Seven combinations of covariates were considered, and their predictive abilities were evaluated by repeated random sub-sampling and herd cross-validations (CV). RESULTS: Addition of the on-farm effects such as herd, days in milk, and parity to spectral data improved predictions as compared to those obtained using the spectra alone. Integrating genomics and/or the top three markers with a large effect further enhanced the predictions. Pedigree data also improved prediction, but to a lesser extent than genomic data. Multiple kernel learning and multilayer BayesB increased predictive performance, whereas PLS did not. Overall, multilayer BayesB provided better predictions than multiple kernel learning, and lower prediction performance was observed in herd CV compared to repeated random sub-sampling CV. CONCLUSIONS: Integration of genomic information with milk FTIR spectral can enhance milk protein trait predictions by 25% and 7% on average for repeated random sub-sampling and herd CV, respectively. Multiple kernel learning and multilayer BayesB outperformed PLS when used to integrate heterogeneous data for phenotypic predictions.

Genome-wide scan highlights the role of candidate genes on phenotypic plasticity for age at first calving in Nellore heifers.

Age at first calving (AFC) plays an important role in the economic efficiency of beef cattle production. This trait can be affected by a combination of genetic and environmental factors, leading to physiological changes in response to heifers' adaptation to a wide range of environments. Genome-wide association studies through the reaction norm model were carried out to identify genomic regions associated with AFC in Nellore heifers, raised under different environmental conditions (EC). The SNP effects for AFC were estimated in three EC levels (Low, Medium, and High, corresponding to average contemporary group effects on yearling body weight equal to 159.40, 228.6 and 297.6 kg, respectively), which unraveled shared and unique genomic regions for AFC in Low, Medium, and High EC levels, that varied according to the genetic correlation between AFC in different EC levels. The significant genomic regions harbored key genes that might play an important biological role in controlling hormone signaling and metabolism. Shared genomic regions among EC levels were identified on BTA 2 and 14, harboring candidate genes associated with energy metabolism (IGFBP2, IGFBP5, SHOX, SMARCAL1, LYN, RPS20, MOS, PLAG1, CHCD7, and SDR16C6). Gene set enrichment analyses identified important biological functions related to growth, hormone levels affecting female fertility, physiological processes involved in female pregnancy, gamete generation, ovulation cycle, and age at puberty. The genomic regions highlighted differences in the physiological processes linked to AFC in different EC levels and metabolic processes that support complex interactions between the gonadotropic axes and sexual precocity in Nellore heifers.

Multitrait meta-analysis identified genomic regions associated with sexual precocity in tropical beef cattle.

Multitrait meta-analyses are a strategy to produce more accurate genome-wide association studies, especially for complex phenotypes. We carried out a meta-analysis study for traits related to sexual precocity in tropical beef cattle (Nellore and Brahman) aiming to identify important genomic regions affecting these traits. The traits included in the analyses were age at first calving (AFC), early pregnancy (EP), age at first corpus luteum (AGECL), first postpartum anoestrus interval (PPAI), and scrotal circumference (SC). The traits AFC, EP, and SCN were measured in Nellore cattle, while AGECL, PPAI, and SCB were measured in Brahman cattle. Meta-analysis resulted in 108 significant single-nucleotide polymorphisms (SNPs), at an empirical threshold P-value of 1.39 × 10-5 (false discovery rate [FDR] < 0.05). Within 0.5 Mb of the significant SNP, candidate genes were annotated and analyzed for functional enrichment. Most of the closest genes to the SNP with higher significance in each chromosome have been associated with important roles in reproductive function. They are TSC22D2, KLF7, ARHGAP29, 7SK, MAP3K5, TLE3, WDR5, TAF3, TMEM68, PPP1R15B, NR2F2, GALR1, SUFU, and KCNU1. We did not observe any significant SNP in BTA5, BTA12, BTA17, BTA18, BTA19, BTA20, BTA22, BTA23, BTA25, and BTA28. Although the majority of significant SNPs are in BTA14, it was identified significant associations in multiple chromosomes (19 out of 29 autosomes), which is consistent with the postulation that reproductive traits are complex polygenic phenotypes. Five proposed association regions harbor the majority of the significant SNP (76%) and were distributed over four chromosomes (P < 1.39 × 10-5, FDR < 0.05): BTA2 (5.55%) from 95 to 96 Mb, BTA4 (5.55%) from 94.1 to 94.8 Mb, BTA14 (59.26%) from 24 to 25 Mb and 29 to 30 Mb, and BTA21 (5.55%) from 6.7 Mb to 11.4 Mb. These regions harbored key genes related to reproductive function. Moreover, these genes were enriched for functional groups associated with immune response, maternal-fetal tolerance, pregnancy maintenance, embryo development, fertility, and response to stress. Further studies including other breeds and precocity traits could confirm the importance of these regions and identify new candidate regions for sexual precocity in beef cattle.

Genetic correlations between body weight, scrotal circumference and visual evaluation scores in Bos indicus cattle.

The objective was to estimate genetic correlations between body weight (BW), scrotal circumference and visual evaluation scores of body conformation measured at standard ages in Guzerat cattle. All measurements were performed at 205 (weaning age), 365, 450 and 550 days of age; for BW, two additional measurements (at birth and 120 days of age) were realized. The data utilized in this study were retrieved from a database of the Brazilian Association of Zebu Breeders that contained information of registered Guzerat animals born between 1970 and 2013. Genetic parameters were estimated in bi-trait analyses by using Bayesian inference. Genetic correlations between BW at 205 and 450 days of age with other traits were high and positive, whereas the correlations between visual evaluation scores with other traits were moderate. Based on correlations herein obtained, we conclude that selection based on BW results in increased visual scores and scrotal circumference, leading to improvements in productive performance and animals with best body conformation.